Three-Sphere Low Reynolds Number Swimmer with a Cargo Container 
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A recently introduced model for an autonomous swimmer at low Reynolds number that is com- 
prised of three spheres connected by two arms is considered when one of the spheres has a large 
radius. The Stokes hydrodynamic flow associated with the swimming strokes and net motion of this 
system can be studied analytically using the Stokes Green's function of a point force in front of a 
sphere of arbitrary radius R provided by Oseen. The swimming velocity is calculated, and shown 
to scale as with the radius of the sphere. 
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Since the pioneering work of G. I. Taylor [l[ we know 
that swimming at low Reynolds number is a nontriv- 
ial task, which is also echoed in the incredibly sophisti- 
cated mechanisms involved in swimming of bacteria 0, Q . 
While nature could overcome the difficulties of swimming 
at low Reynolds number by taking advantage of a contin- 
uum of degrees of freedom, this will be much more chal- 
lenging for an artificially engineered model, as it could 
often use a finite number of degrees of freedom like most 
manmade devices. The relation between the number of 
degrees of freedom and swimming at low Reynolds num- 
ber was considered by Purcell in his seminal work [4|, 
who concluded that one degree of freedom is not enough 
and the minimum required is two. He further elaborated 
on this by proposing a simple swimmer model made of 
three joined rods, whose motion was fully analyzed only 
recently [H due to the complicated nature of the hydro- 
dynamic problem of three rotating finite rods 6] . 

If we take advantage of the positional degrees of free- 
dom rather than the orientational ones, we can construct 
swimmer models that can be much more easily analyzed. 
This is the principle used in the design of a simple swim- 
mer that is made of three spheres connected by two arms 
that open and close Q • The simplicity of the three-sphere 
swimmer model allows us to study various characteristic 
properties of the swimmer itself — such as the effects of 
external load and noise, stress distribution, power con- 
sumption, etc. — in great details [1, Moreover, it also 
makes it possible to study the interaction between two 
[l^ or more of these swimmers, which could be used 
towards a systematic study of the floe behavior of a sus- 
pension of such swimmers [ll|) El, [3 [13] ■ ^ similar 
principle has been used to propose other swimmer mod- 
els that can also be easily analyzed [l^, [l^, [l^. This 
class of swimmer models uses a finite number of compact 
degrees of freedom, and is to be contrasted from those 
that use a continuum [H, [l^ . An artificial microscopic 
swimmer has been recently made using magnetic beads 
connected by DNA linkers, and has been shown to be able 
to propel itself via beating strokes that are induced by an 



oscillating magnetic field [20] ■ It is also possible to design 
swimmers with no moving parts by taking advantage of 
phoretic phenomena, as has been proposed and studied 
theoretically [2l|, [H, [2^ and also realized experimentally 
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The analysis of the three-sphere swimmer model has so 
far been made in the limit where the radii of the spheres 
are the smallest length scales in the system, namely, they 
are much smaller than the distances between the spheres. 
In practical applications, however, one can imagine that 
the swimmer could be used to carry a cargo, which would 
be possible if one of the spheres was blown into a large 
hollow spherical shell, with a radius that is compara- 
ble to the distance between the spheres. A similar ge- 
ometry could also be useful for experimental testing of 
the swimmer design, as a large spherical bead as a part 
of the swimmer would allow optical probing of its pro- 
pelled motion. With these motivations, here we study 
the swimming of the three-sphere model with one of the 
spheres having a finite radius comparable to the other 
length scales. It is possible to fully analyze the motion of 
this model swimmer along the same lines as the simpler 
version, thanks to the Stokes Green's function of a point 
force in front of a sphere of arbitrary radius calculated 
by Oseen in 1927 using the method of images lil]. We 
find that the average swimming velocity scales as 
with the radius of the container. Our results show that 
a simplistic view of "estimating" the velocity of an au- 
tonomous swimmer by dividing the propelling force in 
the system by its Stokes friction coefiicient is incorrect, 
and a only proper analysis of the hydrodynamic flow due 
to the deformations of the swimmer should be used to 
determine the propulsion velocity. 

Consider the geometry of Fig. [T] where one of the 
spheres in the swimmer is in the form of a large con- 
tainer of radius R and the other two have radii ai and 
02. To calculate the swimming velocity of the machine, 
we work in the frame of reference that is co-moving with 
the large sphere. Assuming that there are two force cen- 
ters located at xi and X2 exerting the forces fi and £2, 
respectively, we can write the velocity of the fluid at any 
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FIG. 1: Schematics of the swimmer with a cargo container. 
The instantaneous distance between the two small spheres 
and the center of the large sphere are Ri{t) and R2{t), and 
there three spheres are assumed to always be positioned along 
the z-axis. 
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where V is the velocity of the large sphere in the labora- 
tory frame, and Gij(x, x') is the Oseen Green's function 
that vanishes on the surface of the large sphere [13, . 

Assuming that the two small spheres are located on the 
X-axis (see Fig. [T| at xi = (0, 0, -Ri) and X2 — (0, 0, R2), 
symmetry requires that the device swims along the z- 
direction, namely, V = (0,0,1^). Using no-slip bound- 
ary condition on the two small spheres, and defining the 
quantities 
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we can find their velocities as 
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where Hap = Gzz{^a,^0) for a,f3 = 1,2, and sum- 
mation over repeated indices is understood. Note that 
H12 = 7^21 due to the symmetry of the Green's func- 
tion. Assuming that the small spheres have prescribed 
velocities = Ra, we can determine the forces as 



fa - Sttt? H-'p [Rp + CpV), 



(4) 



where 7Y~^ represent the components of the inverse ma- 
trix for Tiap, with respect to the a,p indices. Note that 
strictly speaking Faxen's theorem [23| requires that an 
additional contribution proportional to the Laplacian of 
the velocity field is added to the expression for the ve- 
locity of the small spheres. However, these contributions 
will generate terms that are higher order in the ratio be- 
tween the sizes of the small spheres and the distances in 
the system, and can hence be neglected here. 

We can now impose the condition of zero external force 
or the vanishing of the total Stokeslet strength as [s^l 



Cafa + Qt^vRV = 0, 



(5) 



which yields the final equation. The above equation can 
be understood as the vanishing of the sum of the forces 
exerted on the fluid at the outer surface of the sphere. 
When we have extended objects rather than just point 
forces, we should enforce the force-free condition by look- 
ing at the stress field in the vicinity of the swimmer and 
show that its integral over any closed surface vanishes. 
In our case here, this stress will have a contribution from 
the sphere, which yields the familiar Stokes force term in 
Eq. (O, plus contributions that have propagated from 
the point forces that are located at distances i?i and i?2. 
This form of the force-free constraint can be related to 
the case of three point-like spheres, where the force-free 
relation reads /i + /2 + /a = 0. To this end, we can sim- 
ply take the limit of i? ^ in Eq. ([5]) and note that the 
viscous force term will produce in this limit as V is the 
velocity of the third sphere. Equation ([5]) proves that it 
is incorrect to assume that the velocity of the bead can be 
extracted from a force balance between the contributions 
from the two force centers and the friction force at the 
bead, which would have (falsely) implied the vanishing 
of Ea + 67rr?i?y. 

Putting in the forces from Eq. ([3]), Eq. ([5]) can be 
solved to yield to propulsion velocity as 



y = -- 



CaTi^pRp 
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(6) 



We also feed this result back into Eq. @ to find the 
forces at the two small spheres. We find 



fa = Snij 



'H^lRp 



(7) 



where "f,S — 1, 2. 

To calculate the swimming velocity we need the ex- 
plicit expression for Gzz(z, z'). which is presented in the 
1927 monograph by Oseen [23|. It reads 



Gz4z,z') = 
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which yields the following expressions for the components 
of the ?i-matrix 
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where the terms involving oi and 02 appear with a dif- 
ferent prefactor according to standard treatments of the 
Oseen tensor to extract the friction coefficient. Putting 
the explicit forms of the TLap matrix elements as func- 
tions of the instantaneous separations Ri (t) and R2 {t) in 
Eq. ([S]), we find a closed form expression for the instan- 
taneous swimming velocity. 

It is instructive to study the asymptotic limits of Eq. 
([5]) in various limits. In the limit where (ai,a2,R) <C 
{Ri, R2, Ri — R2), we can use i?2 = ^2 + U2 and Ri = 
(1 + £2 + ui + U2 and expand the expression in powers of 
Ua- Keeping only terms up to the leading order, we find 
the average velocity as 
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where the averaging over a complete swimming cycle with 
period T is defined by y = ^ dtV{t). This is identical 
to the result obtained previously [1] using the standard 
form of the Oseen tensor, which is a good check for our 
result. (The apparent sign difference with the previous 
result [3| is due to a different convention for the direc- 
tion of the z-axis.) Note that in this case the swimming 
velocity scales as the inverse of the second power of the 
largest length scale in the system, namely 

In the hmit where (ai, 02) ^ (i?i — R, R2 ~ R) ^ R, 
we can use Ri = R + Li + L2 and R2 ~ R + L2 (see Fig. 
[T]), and expand the expression in powers of La/R- We 
find the instantaneous velocity as 
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to the leading order. Expanding for small deformations 
using Li — £1 + ui and L2 = £2 +U2, yields 

- 9 /aia2\ £l (34 + 2^2) 



4 V i?3 J ^2 (^^ 
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{uiU2 - U2U1). (14) 



Equations ([13]) and (fT4|) are proportional to 1/R^, which 
is in marked contrast to the 1 scaling in the case of a 
(more) symmetric swimmer. The difference in the scaling 
suggests that it might be more efficient to distribute the 
size of the swimmer across the three spheres, rather than 
concentrate all of it on one sphere. Such a statement, 
however, needs to be checked in the limit of a symmetric 
swimmer with spheres of large radii compared to their 
surface-to-surface separations. Equation (fT4|) also shows 
that it is beneficial for the two small spheres to be close to 
each other and relatively farther away from the surface of 
the sphere, i.e. £2 ^ £i- Note, however, that this result 
[Eq. (Ull)] is only vahd for (ai, 02) < (^i, ^2) < -R, which 
means that the previous prescription for creating optimal 
swimming efficiency should be written as (ai , 02) ^ ^1 ^ 
£2 ^ R rigorously speaking. 



The present calculation can be readily generalized to 
the case of one large sphere with many small spheres. 
For a line-up of all the spheres along the z-axis, Eqs. 
@, ®, O, ®, ®, and dH) wiU hold true provided 
the indices are run over the number of small spheres in 
the system. For more complicated assortments of small 
spheres, one can still use this formalism provided the 
right tensorial structure is maintained along the way. A 
useful particular case could be a line-up of small spheres 
that oscillate laterally, i.e. perpendicular to the z-axis, 
which would mimic the beating motion of a (finite) tail. 
It will also be interesting to consider the possibility of 
a number of such structures attached to the bead inter- 
acting with each other hydrodynamically. These "active 
tails" could potentially synchronize with each other un- 
der certain conditions, to generate an efficient propul- 
sion. Our systematic approach might be able to provide 
useful insight into the intriguing phenomenon of collec- 
tive hydrodynamic patterns and phases [3l|, [H, [s^ . One 
can also study the efficiency of the swimmer in terms of 
power consumption and useful work, which for this swim- 
mer follows very closely similar studies that are presented 
elsewhere [8j]. 



In conclusion, we have generalized the simple model 
of a low Reynolds number swimmer to the case where 
one of the spheres is large, and analyzed its motion using 
the Oseen tensor of a point force in front of a sphere 
which imposes no-slip boundary condition on the Stokes 
flow. The results found here could help us in possible 
experimental probes of self-propelled devices. 
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